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The available results from the inelastic neutron scattering experiment performed on the quasi-two 
dimensional spin | anti-ferromagnetic material La 2 Cu 04 have been analysed theoretically. The for¬ 
malism of ours is based on a semi-classical like treatment involving a model of an ideal gas of mobile 
vortices and anti-vortices built on the background of the Neel state, using the bipartite classical spin 
configuration corresponding to an XY- anisotropic Heisenberg anti-ferromagnet on a square lattice. 

The results for the integrated intensities for our spin | model corresponding to different tempera¬ 
tures, show occurrence of vigorous unphysical oscillations, when convoluted with a realistic spectral 
window function. These results indicate failure of the conventional semi-classical theoretical model 
of ideal vortex/anti-vortex gas arising in the Berezinskii-Kosterlitz-Thouless theory for the low spin 
magnetic systems. A full fledged quantum mechanical formalism and calculations seem crucial for 
the understanding of topological excitations in such low spin systems. Furthermore, a severe dis¬ 
agreement is found to occur at finite values of energy transfer between the integrated intensities 
obtained theoretically from the conventional formalism and those obtained experimentally. This 
further suggests strongly that the full quantum treatment should also incorporate the interaction 
between the fragile-magnons and the topological excitations. This is quite plausible in view of the 
recent work establishing such a process in XXZ quantum ferromagnet on 2D lattice. The high spin 
XXZ quasi-two dimensional antiferromagnet like MnPSs however follows the conventional theory 
quite well. 

Keywords: Spin dynamics. Topological spin excitations, Berezinskii-Kosterlitz-Thouless transition, 

Spin 1/2 easy plane Anti-ferromagnet. 

Highlights: 

• Inadequacies in the conventional meron gas phenomenology in explaining the spin dynamics 
corresponding to layered low-spin anti-ferromagnets. 

• Requirement of a full fledged quantum mechanical formalism for the understanding of spin 
dynamics induced by both topological and conventional excitations. 

• Necessity of a proper understanding of the interaction between the conventional and 
topological excitations at the quantum level. 


I. INTRODUCTION 

Spin dynamics in low dimensional magnetic systems 
have generated a significant research interest during the 
last three decades [TH5T]. Different types of one and quasi- 
one dimensional as well as two and quasi-two dimensional 
systems have been studied both experimentally and the¬ 
oretically to probe and understand the spin dynamics 
arising from the conventional spin wave excitations and 
their mutual interactions, as well as from vortex/meron 
and soliton like topological excitations [l5HT7llT^[2(711?7] . 
In many of these systems the topological excitations of 
soliton and vortex/meron type occur naturally as they 
are thermodynamically feasible. 

Motivated by the distinct possibilities of applica¬ 
tions towards building of magnetic devices, the quasi-two 
dimensional systems have attracted a renewed research 
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interest in recent times. Magnetic vortices present in 
these systems, have proved to be potential candidates for 
switching devices [821 - 135] . Direct experimental evidences 
of such vortices have been verified by both the Magnetic 
Force Microscopy (MFM) and the spin-polarized Scan¬ 
ning Tunnelling Microscopy (STM) [8611^ . 

The spin dynamics in many of the above magnetic 
systems has been investigated experimentally using the 
Inelastic Neutron Scattering (INS) and the Nuclear Mag¬ 
netic Resonance (NMR) techniques. These include quasi- 
one dimensional systems such as CsNiF^, layered sys¬ 
tems such as K 2 CUF 4 , Rb 2 CrCl 4 , LiCr 02 , magnetically 
intercalated graphites, layered ruthenates, layered man- 
ganites and the high-Tc cuprates [TBHT^ |3I] ■ In 

the INS experiments performed on several of above ma¬ 
terials, the existence of a prominent “central peak” (at 
hoj = 0) has been confirmed in the plot for the dynami¬ 
cal structure function S'(q, uj) vs. neutron energy transfer 
‘huj’ in the constant ‘q’ scan [2111^. These findings fur¬ 
ther serve as the motivation behind the huge variety of 
experiments performed on the layered magnetic systems. 
Moreover, the advancement of numerical and computa- 
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tional techniques also contributed to the understanding 
of the possible role of both the spin waves and the topo¬ 
logical excitations in emergence of the central peak. 

Kosterlitz and Thouless, and Berezniskii indepen¬ 
dently introduced the concept of vortex and anti-vortex 
like topological spin excitations in the two dimensional 
classical magnetic (spin) systems [55MD] . According to 
their ideas, there exists a non-conventional topological 
phase transition (known as Berezinskii- Kosterlitz- Thou¬ 
less or BKT transition) characterized by the crossover be¬ 
tween binding to unbinding phases of vortex- anti-vortex 
pairs at a transition temperature Tbkt- Below this tem¬ 
perature all the vortices and anti-vortices are in a bound 
state and above this temperature some of them start 
moving freely. Further analytical and numerical stud¬ 
ies and suitable extension of these approaches led to the 
proposal for the existence of topological vortices and anti¬ 
vortices in pure XY model and merons and anti-merons 
in XY- anisotropic Heisenberg model, for both ferro¬ 
magnetic and anti-ferromagnetic types, on two- dimen¬ 
sional lattices |1TI - HS] . Furthermore, approximate ana¬ 
lytical calculations and Monte Carlo-Molecular Dynam¬ 
ics (MCMD) simulations have strongly suggested that 
the freely moving topological excitations in the regime 
T > Tbkt, contribute non-trivially to the spin-spin cor¬ 
relation and give rise to the “central peak”, as mentioned 
above [¥TI - Hi] . The occurrence of such a “central peak” 
in quasi-two dimensional magnetic systems is now unan¬ 
imously believed to be the signature for the dynamics of 
mobile topological excitations in a layer. 

In spite of a lot of studies on the anisotropic Heisen¬ 
berg model on two dimensional lattices the contributions 
of vortices/anti-vortices to the dynamics of the model, 
especially the detailed quantitative features of the DSFs 
are still not completely understood. 

In the context of two-dimensional magnetism, the 
undoped (anti-ferromagnetic and non-superconducting) 
phases of the high Tc cuprate systems are believed to 
be excellent examples of two dimensional XY anisotropic 
Heisenberg Hamiltonian in an appropriate temperature 
regime. One member of this class of systems is La 2 CuOi, 
on which extensive INS experiments have been per¬ 
formed 1301 [31]. This is a truly spin-1/2 layered anti- 
ferromagnet. The intra-layer integrated intensity corre¬ 
sponding to the results of INS experiment performed on 
La 2 Cu 04 , exhibits a central peak when plotted against 
the neutron energy transfer Huj (or frequency ‘w’). 

It has been shown that the results of vortex gas 
phenomenology and numerical simulations lead to an 
anomaly in the case of layered anti-ferromagnetic sys¬ 
tems having very low spin values (S=l/2) [37|. Strik¬ 
ingly enough, the value of Tbkt obtained from Renor¬ 
malization group analysis and numerical simulations is 
four (4) times the value of Tbkt calculated from the 
classical expression obtained by Kosterlitz and Thouless 
[47j . Furthermore, in our previous work we have already 
established that for quasi-two dimensional ferromagnetic 
systems having low spin values (S= 1/2) the conventional 


semi-classical like treatment involving the ideal gas of 
unbound vortices/merons and anti-vortices/anti-merons 
corresponding to high temperature regime T > Tbkt, 
shows large inconsistency with the experimental situation 
and exhibits unphysical behaviour |1S|. In this case the 
theoretical dynamical structure function (DSF) turns out 
to be negative for a wide range of energy transfers! How¬ 
ever the range, over which the theoretical DSF remains 
positive, increases when the value of the spin is increased 

[48] . These facts motivate us to investigate and test in de¬ 
tail the applicability of the semi-classical-like treatment 
mentioned above to the INS results corresponding to real 
anti-ferromagnetic systems with S=l/2. For this exercise 
we select La 2 Cu 04 as the reference system [3(71 l3T] . 

It is worthwhile to mention that the BKT transi¬ 
tion can also be identified at the transition tempera¬ 
ture Tbkt where the spin-stiffness jumps discontinuously 
from a universal value below Tbkt to zero above Tbkt 

[49] . Moreover, in anti-ferromagnetic model possessing 
Ising like anisotropy (in the z-direction) on two dimen¬ 
sional lattices, with external field being applied in the z- 
direction, such a discontinuous jump has been observed 

[50] . However, the discontinuous jump in this case may 
have its origin different from the vortex-anti-vortex un¬ 
binding mechanism since it is well known that the BKT 
transition in magnetic systems can only occur when the 
anisotropy is XY like [411 - 144] . 

The plan of the paper is as follows:- in section H we 
describe the formulation of our semi-classical like treat¬ 
ment; in section HI we discuss our calculations and re¬ 
sults and finally in section IV the conclusions and discus¬ 
sions of our present investigation are presented. 


II. MATHEMATICAL FORMULATION 

The dynamics of mobile topological excitations in an 
anti-ferromagnetic system on a two-dimensional square 
lattice have been analysed both analytically and numeri¬ 
cally [3T1H5] . The analytical studies have been performed 
by assuming a classical ideal gas of vortices/merons 
where the vortices/merons obey Maxwell’s velocity dis¬ 
tribution. The model system is described the XY- 
anisotropic Heisenberg (XXZ) Hamiltonian, viz., 

n = -j + sf^sy^ + xs:^s;^), (i) 

{ihpq) 

where {ij,pq) label the nearest neighbour sites on 
a two-dimensional square lattice and J(< 0) is the 
anti-ferromagnetic exchange coupling. Here A is the 
anisotropy parameter whose pure XY and isotropic 
Heisenberg limit correspond to A = 0 and 1 respectively. 

The structures of the vortices/merons have been ob¬ 
tained by solving the classical equations of motion corre¬ 
sponding to the Hamiltonian given by eqn. 0. In deriv¬ 
ing the classical equations of motion the spins have been 
considered to be classical objects (classical spin fields 
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S{r,t)) as a function of position coordinates and time, 
which are defined on the entire lattice. At even or odd 


lattice sites these spin fields become identical to the fol¬ 
lowing bi-partite spin configurations: 


gepn ^ + Oij) cos($ij + sin(0ij + 0^^) sin($y + cos(0ij + %))], 

S°f^ = -S[sm{Qij - Oij) cos($ij - (t>ij), sin(0ij - 0^) sin($y - ), cos(0ij - %))], (2) 


where ’’even' and ’’odd’ signifies the two different sub¬ 
lattices m- The static spin configuration corresponding 
to the merons are described by the capital angles 0(r) 
(polar) and <i)(r) (azimuthal), and the time dependent 
small angles 0{r,t) and (/>(r,<) describes the correspond¬ 
ing deviations from the static structure due to the motion 
of the merons and the spin dynamics above BKT transi¬ 
tion temperature [HI 03] . The expression of the vortex 
core radius is given by [H |45] 


From the above considerations the in-plane dynam¬ 
ical structure function (in-plane DSF) w) is given 

by, 


5""(q,w) 


>g(g + l) 

27r (a; 2 + 72 [i + (^q*)2])2’ 


(4) 


with 7 = where q* = (qo —q); qo = (tt/o, tt/o) and 
in our case S' = 4. The above expression for the in-plane 
dynamical structure function is a squared Lorentzian ex¬ 
hibiting a central peak at a; = 0 in ‘w’-space for constant 
q- scan and exhibiting a central peak at the zone bound¬ 
ary of the first Brillouin Zone (BZ) in the ‘q’ space for 
constant w-scan [HI |3S]- In the above expression u is 
the root mean square (rms) velocity of the vortices and 
is given by, 

u = (5) 

where nl ~ (2^)“^ is the density of free vortices at T > 
Tbkt mi - Here ^ is the intra-layer two-spin 

correlation length due to the presence of vortices, where 
^0 is of the order of lattice spacing ‘o’; r = — 

1) is the reduced temperature and 6 is a dimensionless 
parameter whose numerical value is generally around 1.5 
[221 Si] ■ The quantity S"^’'(q, w) is sensitive to the in¬ 
plane structure of the vortices/merons [^ I45j . 

Again the effective analytical expression for the out- 
of-plane dynamical structure function (out-of-plane DSF) 
S'^^(q, w) in the limit of very small ‘q’ is given by pUfiS] 




32(1+ A)V2 ^uq 


( 6 ) 


at a; = 0, when plotted in the constant q-scan. The 
function 5'^^(q, w) is sensitive to the out-of-plane shape 
of the vortices/merons [HUS]. 

In the case of layered systems, in a suitable regime 
in the parameter space comprising of temperature and 
wave vector where these systems behave effectively as 
two-dimensional systems, the integrated intensity corre¬ 
sponding to a typical inelastic neutron scattering exper¬ 
iment is given by 

4(^) = ['^s:°‘°'{^2r>,i^)dqxdqy, (7) 

a 

where the quantity 5'““(q2D,w) represents the intra¬ 
layer in-plane spin dynamical structure function when 
a = X and y and the intra-layer out-of-plane spin dy¬ 
namical structure function when a = z |1511521153j . 

The approach in our present work is quite similar to 
that adopted earlier in the case of ferromagnetic systems 
[48] . In order to compare theory with experiment the dy¬ 
namical structure function, obtained from the model un¬ 
der consideration, is multiplied with the resolution func¬ 
tion R{t) (in time domain) or convoluted with R{uj — uj') 
(in the frequency domain), as has been done in the case 
of ferromagnetic system [48l [52] (54] . Hence, the compo¬ 
nents of the convoluted integrated intensity comes out to 
be. 


TOiOi 

-^conv 


(w) 


/ 


dqx dqy 


dw'i?(w-a;')S'““(q2D,w'). (8) 


The resolution function has to be chosen in such a way 
that minimum ripples occur at the end points of the res¬ 
olution width. The different parameters of the resolu¬ 
tion function can be obtained from the resolution half 
width or the full width at the half maximum (FWHM) 
which are quoted in the experiments. Since X and Y 
components of the spins are symmetric i.e., S'’'’^(q, w) = 
^^^(q, uj) the total intensity comes out to be, 

I{uj) = 2r^{oj) + P^{oj). (9) 


Furthermore, keeping in mind the low spin situa¬ 
tion the quantum mechanical detailed balance condition 
is incorporated in our formalism [55] . Then semi-classical 
estimate for I{uj), denoted by Isci^^) is recovered by the 
relation. 


2 

l + exp{f^) 


4(w), 


The above form of the out-of-plane dynamical structure 
function is a Gaussian, exhibiting again a central peak 




( 10 ) 
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where the factor ^ is the detailed balance factor 

l+eap(j^) 

and is called the Windsor factor [56l [57] . The superscript 
‘SC’ stands for the term semi-classical. 

As has been pointed out in the case of ferromagnetic 
systems, this approach of ours is truly ‘semi-classical-like’ 
in the sense that the expression for the rms vortex veloc¬ 
ity ‘u’, as given in eqn. (4), contains ‘/i’ and in addition 
the quantum mechanical detailed balance condition has 
been incorporated through the Windsor factor [48], 

It is worthwhile to mention that the above formu¬ 
lation based on dilute vortex/meron gas phenomenol¬ 
ogy hold for unbound anti-vortices/anti-merons too on 
the basis of the assumption that the vortices/merons 
and anti-vortices/anti-merons do not interact with each 
other. 


III. CALCULATIONS AND RESULTS 


In this work the formalism of Section II is applied to an 
anti-ferromagnetic material La 2 CuOi on which inelastic 
neutron scattering experiments (INS) involving polarized 
neutron beam have been performed [sniEi]. The ma¬ 
terial is an XY-anisotropic quasi-two-dimensional spin 
1/2 quantum Heisenberg anti-ferromagnet. The mag¬ 
netic lattice structure of it is composed of stacking of 
two-dimensional square lattices [sniisi]. The spin Hamil¬ 
tonian relevant to the above material is given by, 


n = {-j s* ■ s, + JaJ2 5 ] s, • s, 


<i,3> 


<i,3> 


<i,k> 


intra-layer part 


inter-layer part 

( 11 ) 

where < i,j > represents the intra-layer nearest neigh¬ 
bour interaction and < i,k > represents the inter-layer 
nearest neighbour interaction. In the above Hamiltonian, 
J is the isotropic part and Ja is anisotropic part of the 
intra-layer exchange coupling and, J' is the inter-layer 
exchange coupling. The ordering temperature i.e., the 
Neel temperature for the quasi-two dimensional system 
{La 2 Cu 04 ) is given by Tm = 240 K. The intra-layer part 
of the above Hamiltonian (111 can be simplified, by ex¬ 


pressing {J — Ja) as AJ, to obtain the model Hamiltonian 
0. where A is the anisotropy parameter. The relevant 
physical parameters corresponding to La 2 CuOi are given 
in the TABLE |!] [58]. 


Parameter 

Magnitude 

J (intra-layer) 

~ 1345 K 

Ja (intra-layer) 

~ 0.269 K 

J' (inter-layer) 

~ 0.04 K 

anisotropy parameter (A) 

0.9998 

lattice parameter(a) 

5.39 A 

Neel temperature (Tjv) 

240 K 


TABLE I. Relevant parameters for La 2 Cu 04 


Next we try to determine the temperature range over 
which the material La 2 CuOi behaves effectively as a two- 
dimensional material. From the neutron scattering data 
for quasi-two-dimensional spin 1/2 XY-anisotropic fer- 
romagnet K 2 CuFi, it was found that in a temperature 
regime Ti < T < T 2 , where the lower (Ti) and the upper 
(T 2 ) limits are defined by the following relations, 

m) = 



the system behaves as a 2D XY-like anisotropic system 
[mill]. Assuming that the above phenomenological ar¬ 
gument holds for the layered anti-ferromagnetic systems 



DataFR 



Temperature T (In K} 


FIG. 1. Fitting of the experimentally obtained inverse correla¬ 
tion length with the corresponding theoretical expression (see 
eqn. {1^ ). Solid line corresponds to the theoretical expres¬ 
sion. BKT transition temperature is Tbkt = 270A. 

as well, we determine the above two temperature limits 
as Ti « 260A and T 2 « 360A for La 2 CuOi. Within 
this temperature regime the Hamiltonian ( jllj) can ef¬ 
fectively be represented by the Hamiltonian (|H. It is 
worthwhile to point out that since in the above temper¬ 
ature regime the system is effectively a two-dimensional 
one, long range anti-ferromagnetic ordering is absent in 
this regime |59j . Further, within the above mentioned 
temperature range the BKT inspired ideal vortex/meron- 
gas phenomenology is valid and we can therefore use the 
theoretical expression for the inverse correlation length 
(expressed in r.l.u), 

k{T) = (13) 

TT 

as predicted by Kosterlitz and Thouless, to fit the exper¬ 
imentally obtained inverse correlation length [38] . This 
gives the value of the BKT transition temperature as 
Tbkt ~ 270A for La 2 CuO^ (see FIG. [^. In this work 
we shall make use of the above value of Tbkt to calculate 
the convoluted in-plane integrated intensity, (w) and 
the convoluted out-of-plane integrated intensity 
by making use of eqns. (|^ to ([^, and then eqns. 0 and 

0 . 
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We now study the convoluted in-plane integrated in¬ 
tensities Iconvi^) different temperatures. The expres¬ 
sion for the Iconvi^) is given by eqn. (§ with a = X, 
where the in-plane DSF, S'“'^(q 2 D,w) is given by eqn. 
Q. In the experimental investigations on La 2 CuOi, to 
find the neutron intensity as a function of momentum 
transfer ‘fi.q’ the scans in the ‘q’-space have been per¬ 
formed about the zone boundary of the first BZ within 
the range, —0.1 < q* < 0.1, expressed in r.l.u [SnilSI]- In 
calculating using eqn. (13) we have also made 

use of the above mentioned regime only. The resolu¬ 
tion function has been chosen in the form of the Tukey 
window to convolute the in-plane DSF. This is one of 
the most commonly used spectral smoothing functions 
in the field of spectral analysis [SOI m]. The experi¬ 
mental resolution width is 1.4 meV at the full width at 
half maximum (FWHM) as specified in the experiment 
[30l We compute for four different tem¬ 

peratures, viz., 290 K, 320 K, 350 K and 375 K. The 


unwanted oscillations diminish considerably when we use 
the modified Tukey function. More interestingly, FIG. 




FIG. 2. The plot of convoluted (with Tukey window func¬ 
tion) in-plane integrated intensity Iff„^{uj)\sc at four dif¬ 
ferent temperatures, viz., 290 K, 320 K, 350 K and 375 
K. The rms velocities at these temperatures are u = 
0.00365^, 0.0836^, 0.085^, 0.2323^ respectively. 


semi-classical convoluted in-plane integrated intensities 
denoted by, sc are plotted as functions of energy 

transfers in FIG. where eqns. © and ( |I0| have been 
used. The figures clearly exhibit that 7c(mjjT^)ls'C' oscil¬ 
lates vigorously after convoluting with the Tukey func¬ 
tion, although the ‘central peak’ still persists. This is 
quite contrary to what we experienced in the case of fer- 
romagnet [15] . 

To avoid such oscillations we later tried perform¬ 
ing the above calculations using a modified version of 
the Tukey function (see eqns. (A4) and (A5) in the Ap¬ 
pendix). The integrated intensity (at 290 K) correspond- 
ing to this new window function is also plotted in FIG 
along with the same corresponding to the use of Tukey 
function. From this figure it is clearly visible that the 


FIG. 3. The plot of convoluted in-plane integrated intensity 
Iconv{^)\sc at 290 K. The red solid line corresponds to the 
Ifif„v{^)\sc obtained by using the Tukey window function (see 
The black solid line corresponds to the 
obtained by using the modified Tukey window function (see 


indicates that at even higher temperature, viz., at around 
375 K (« 1.^88Tbkt) both the Tukey function and the 
modified Tukey function lead to very similar results. The 
oscillations are totally absent in the theoretical plot of 
^conv{^)\sc VS. energy transfers corresponding to both 
the resolution functions. It is worthwhile to mention that 
the corresponding temperature T=375 K (> T 2 ) falls 
just outside the range Ti < T < T 2 within which the 
BKT phenomenology remains valid. However, the use 
of such a modified Tukey function may wipe out some 
of the genuine and intrinsic fluctuations present in the 
anti-ferromagnetic systems in two dimensions. Hence, 
we make use of the Tukey function only for our purpose 
because it is very frequently used in the field of spectral 
analysis. 

We further notice a slight shift in the position of 
the central peak. This is due to the inclusion of detailed 
balance condition. The important point here is that the 
shift is well within the resolution width 1.4 meV at the 
FWHM, and hence the peak is truly a central peak situ¬ 
ated at huj = 0. 

Exactly the same results hold for /^j'„^(a;)|5C' which 
is obvious from the symmetry argument. The normal¬ 
ization factor required for the quantitative comparison 
between the theoretical and the experimental results is 
estimated from the neutron count corresponding to the 
experimental results for La 2 CuOi. 

We now evaluate the out-of-plane integrated inten¬ 
sity Iconv{^)\sc for the same set of temperatures as have 
been considered earlier for the evaluation of Ifonv{^)\sc 
(see FIG [^. The expression for the /®^„(a;) is given 
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FIG. 4. The plot of convoluted in-plane integrated intensity 
Iconv{^)\sc ot 375 K. The red solid line corresponds to the 
Iconv{^)\sc obtained by using the Tukey window function (see 
The black solid line corresponds to the Iff„^{Ljj)\sc 
obtained by using the modified Tukey window function (see 


by, eqn. ® with a = z, where the out-of-plane DSF, 
is given by eqn. <©■ In this case we find 



FIG. 5. The plot of convoluted out-of-plane integrated inten¬ 
sity I((nv(.^)\sc 0,1 four different temperatures, viz., 290 K, 
320 K, 350 K and 375 K. The rms velocities at these tempera¬ 
tures are u = 0.00365;-^, 0.0836;-^, 0.085^, 0.2323^ 

tnat ’ tnat ’ inat ’ ^nat 

respectively. The order of magnitude of the I((nv{oj)\sc ot the 
above mentioned four different temperatures are 
10~® and 10~^ respectively. 

that the out-of-plane integrated intensity oscillates only 
at lower temperatures near T — Tbkt- It is worth re¬ 
calling here that the magnitude of the out-of-plane in¬ 
tegrated intensity is proportional to the density nl and 
the rms velocity u of the free vortices/merons. Since 


both density nl and the rms velocity u increases with 
increasing temperature the out-of-plane part of the spin- 
spin correlation (see eqns. ® and ([^) acquire domi¬ 
nance (considerable magnitude) only at higher tempera¬ 
tures much above Tbkt- 

Furthermore, the absolute magnitude of the inte¬ 
grated intensity is higher (almost 10^ times 

for the highest temperature considered here) than that of 
Iconv{^)\sc at temperatures above Tbkt- This is so be¬ 
cause, Iconv{^)\sc is proportional to nl and it increases 
with the increasing value of the rms velocity u. The 
typical energy scales involved in the dynamics of mobile 
vortices/merons corresponding to the anti-ferromagnetic 
system La 2 CuOi are such that nl is very small (com¬ 
pared to the case of ferromagnet where the free vortex 
number density turns out to be appreciable [48]). For 
the present case of La 2 CuOi at T = SSOFiT, the numeri¬ 
cal value for the density of free vortices comes out to be 
nl = 1.36 X 10“'* and the same for the rms veloc¬ 
ity comes out to be m = 9163 m/ sec = 0.085^^, where 
inat = (« 5 X 10“*® sec) is the natural time scale for 

the system. In contrast, in the case of ferromagnetic sys¬ 
tem K 2 CUF 4 the value of the density of free vortices was 
found to be nl = 1.009 x 10“® a~^ and that for the rms 
velocity was found to be it = 87.07m/sec = 0.1352^^ 
at T = 6.75iF, where tnat = 6.4 x 10“*® sec is the natural 
time scale [48] . 

Interestingly enough, the unbound merons/vortices 
above Tbkt move much faster (with a rms velocity of 
9163 m/ sec at 350 K) than a typical Copper (Cu) atom 
whose rms velocity (generally considered to be the ther¬ 
mal velocity) is around 370.6 m/ sec at 350 K. 

The integrated intensities computed above corre¬ 
spond to the contributions only from the mobile vor¬ 
tices/merons. The experimental data whereas, contain 
contributions from both the mobile vortices/merons and 
fragile “spin wave like” modes. This spin wave like modes 
are damped and largely decaying above the Neel temper¬ 
ature. The extraction of the mobile vortex contribution 
from the experimental data is crucial for a more accu¬ 
rate comparison between the theoretical results and the 
experimental data and to do this one has to subtract 
the contributions from the above mentioned fragile “spin 
wave like” modes from the experimental data. It is worth 
recalling that in the case of ferromagnetic system, the 
fragile mode contributions have been subtracted by as¬ 
suming the fragile mode contributions above Tc to be the 
same as the usual spin wave contributions below Tc- This 
assumption is however valid if and only if the tempera¬ 
ture under consideration is in the near vicinity of the 
Curie temperature (Tc) of the system |35|. 

In the present case corresponding to the anti¬ 
ferromagnetic system La 2 Cu 04 , the situation is some¬ 
what different from the case of ferromagnetic systems in 
the sense that the temperatures dealt with are far above 
the Neel temperature Tjy. Hence the above mentioned 
procedure, which was followed for the ferromagnetic sys- 
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terns to subtract the fragile mode contributions, is not 
valid. 

Moreover at any finite temperature above Tbkt, it is 
to be kept in mind that not all vortices/merons are freely 
moving and that bound vortex-anti-vortex pair density 
remains finite. Hence, one has to estimate further the 
contribution from these bound vortex-anti-vortex pairs 
at different temperatures and subtract them from the 
experimental data. To estimate the bound vortex contri¬ 
bution we have tried to apply the same methodology that 
has been outlined and used earlier for the ferromagnetic 
system [35]. However, unlike the case of ferromagnet 
in this case the methodology leads to an unphysical be¬ 
haviour viz., the vanishing of the out-of-plane DSF (see 
eqn. §) in the limit m —> 0. Hence the estimation of 
the bound vortex contribution has not been carried out 
here. 

The above calculations for the components of inte¬ 
grated intensities enable us to estimate theoretically the 
total integrated intensity using In FIGs. d and 



FIG. 6. The plot of convoluted total integrated intensity 
Iconv{ia)\sc at 290 K. The red solid line corresponds to the 
Iconv{ia)\sc o btain ed theoretically by using the Tukey window 
function (see ( A4)J- The dots are the experimental data. 


[^the total intensity Iconviy>)\sc has been compared to 
the experimental data at two different temperatures, viz., 
290K and 350K respectively. The contributions from the 
fragile “spin-wave-like” modes can not be filtered out 
for the reasons stated earlier. It is clear from FIG. [6| 
that the total intensity Iconvito)\sc also oscillates vigor¬ 
ously at both the temperatures, when convoluted with 
the Tukey function. Besides, the theoretical results show 
negative values for Iconvi^)\sc with both the resolution 
functions! Moreover, the magnitude of the total intensity 
Iconv{^)\sc^ obtained theoretically at finite energy trans¬ 
fers, is very far from the corresponding values obtained 
in the experiment. 

The inclusion of quantum mechanical detailed bal¬ 
ance factor in our semi-classical like treatment has again 
caused a shift in the position of the central peak of the in¬ 



energy transfer h(o/27t in meV 


FIG. 7. The plot of convoluted total integrated intensity 
Iconv{a})\sc at 350 K. The red solid line corresponds to the 
Iconv{a})\sc obtained theoretically by using the Tukey window 
function (see egn. (Af)). The dots are the experimental data. 


tegrated total intensity at both the temperatures. How¬ 
ever, this shift is well within the resolution width and 
therefore it is a genuine central peak at zero energy trans¬ 
fer. 


It is worthwhile to mention that all the above results 
based on dilute vortex/meron gas phenomenology hold 
for unbound anti-vortices/anti-merons too. 

Let us further calculate the zeroth moment of 
the semi-classical dynamical structure function (DSF) 
Sg°Q'"{q,uj) (or simply the moment) using the following 
formula HD, 



S{S + 1), 

S{S + 1) 

(14) 


where is given by equation (10) and in obtain¬ 

ing the same the integration over the wave vector space 
in equation ([^ is performed over the first Brillouin zone 
(B.Z.) with contributions from both vortices and anti¬ 
vortices being summed; S is the value of the spin corre¬ 
sponding to the system under consideration and in our 
case it is S = 1/2. The above equation signifies that if 
the spin dynamics is entirely captured by the DSF, the 
value of the zeroth moment must be S{S + 1). 















T 

T 

T'bkt 

moment in the 
unit of S{S -1- 1) * 

290 K 

1.074 

0.0400 

320 K 

1.185 

0.2412 

350 K 

1.296 

0.3583 

375 K 

1.388 

0.3770 


TABLE II. The zeroth moment of the semi-classical dynamical 
structure function corresponding to the dynamics of mobile 
vortices and anti-vortices{* corresponding to the use of Tukey 
function.) 

The values of the moment of the semi-classical convo¬ 
luted DSF corresponding to the use of Tukey function at 
four different temperatures are tabulated in TABLE [IT] 
At 290A' i.e. around 1.074:Tbkt, the combined dynamics 
of mobile vortices and anti-vortices capture only about 
4% of the entire spin dynamics of the system. However, 
at higher temperatures around 1.296Tbkt-, the combined 
dynamics of mobile vortices and anti-vortices capture 
more than 35% of the entire spin dynamics of the system . 
This happens because at lower temperatures near Tbkt, 
the number of freely mobile vortices and anti-vortices 
are not large enough to capture the whole spin dynam¬ 
ics and the presence of fragile or damped spin waves (or 
single magnons and multi-magnon like modes) makes im¬ 
portant enough contribution to the spin dynamics. At 
higher temperatures however, more topological excita¬ 
tions become free and drive a large portion of the spin 
dynamics of the system. At this point, it is worth men¬ 
tioning that for quantum ferromagnetic systems on two 
dimensional square lattice it has been shown that the for¬ 
mation of topological excitations of vortex/ meron types 
from the fragile magnons and multi-magnon composites 
is quite plausible [S2]. Moreover, some of the collective 
modes (i.e. magnon and multi-magnon modes) are ex¬ 
pected to stay intact with their damped nature and thus 
can provide a significant contribution to the spin dynam¬ 
ics. In analogy with the three dimensional systems where 
above the Curie temperature (Tc) the magnon-like collec¬ 
tive excitations become fragile and damped, for pure two 
dimensional systems (where = 0) the collective excita¬ 
tions become fragile at any finite temperature [561I62H67] . 
This process is expected to be operative too in the anti¬ 
ferromagnetic systems on pure two-dimensional lattices. 

To summarize, we find vigorous oscillations in the 
convoluted in-plane integrated intensity when the Tukey 
window is used. These oscillations vanish only at higher 
temperatures which are outside the regime of validity of 
the BKT phenomenology. The use of a modified or re¬ 
fined Tukey function substantially removes the unwanted 
oscillations in the convoluted in-plane integrated inten¬ 
sity Iconv{^)\sc- Strikingly enough, at T = 350 K (1.296 
Tbkt), we still find negative values of Iconv{^)\sc even 
using the modified Tukey window function. However, 
outside the temperature regime where the BKT phe¬ 
nomenology is valid, computations with both the window 
functions give very similar results for //)f„.j,(w)|sc>- The 


possible explanation for this is that at higher temper¬ 
atures the quantum effects are less prominent even for 
S' = i anti-ferromagnet. Therefore, the modified Tukey 
window may actually be suppressing quantum fluctua¬ 
tion as well quite efficiently. The out-of-plane integrated 
intensities Iconvi^)\sc (computed at different tempera¬ 
tures) are found to be sensitive to the choice of window 
function only at temperatures which are not very far from 
Tbkt (around 1.074 Tbkt)- Furthermore, it contributes 
negligibly to the total integrated intensity Iconv{to)\sc 
and hence the nature of the convoluted total integrated 
intensities at different temperatures turns out to be quite 
similar to that of the convoluted in-plane integrated in¬ 
tensities. However, the detailed quantitative compari¬ 
son between the theoretical results and the experimental 
results corresponding to the total integrated intensities 
Iconv{^^)\sc reveals that even though our “semi-classical 
like” theory is able to predict the occurrence of the cen¬ 
tral peak, the magnitudes of the /conu(w)|sc for finite 
values of energy transfer, obtained from theoretical anal¬ 
ysis, differ by a huge factor from the corresponding exper¬ 
imental values. Moreover, apart from the spin dynamics 
induced by the mobile topological excitations, the frag¬ 
ile magnons and multi-magnon modes are quite likely to 
make important contribution towards this dynamics. 


IV. CONCLUSIONS AND DISCUSSIONS 

It is a well known fundamental fact that the spin- 
dynamical structure function (DSF) and hence the 
integrated intensity corresponding to all real magnetic 
systems must be positive definite [m [521 [53]. Our 
detailed calculations and analysis however, brings 
out an important fact that for quasi-two dimensional 
low-spin anti-ferromagnetic systems, the semi-classical 
treatment of ideal gas of mobile vortices/merons (anti- 
vortices/anti-merons) leads to negative values of total 
integrated intensities even in the low energy regime, 
when convoluted with commonly used resolution func¬ 
tions. Similar results were obtained for 5 = 1/2 layered 
ferromagnets |3H]- Such a behaviour is purely unphys¬ 
ical. These facts highlight the inapplicability of the 
semi-classical vortex/meron gas phenomenology to the 
quasi-two dimensional (or rather layered) low-spin anti¬ 
ferromagnetic systems very strongly and necessitates the 
demand for a full quantum treatment of the problem. 
Let us further point out that the value of huj at which 
the onset of such unphysical behaviour occurs depends 
on the value of the spin (5). It has been shown in the 
case of ferromagnetic systems that the regime over which 
this unphysical behaviour persists, shrinks as the value 
of S increases |48] . Similar behaviour is expected in the 
anti-ferromagnetic systems also. This may be verihed 
from the studies on the compounds based on Gd and Mn. 

To the best of our knowledge there are no mate¬ 
rials based on (spin |) which can be modelled 
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by 2D XXZ Hamiltonian (XY anisotropic Heisenberg 
Hamiltonian). However, some high spin magnetic sys¬ 
tems such as Rb 2 MnF 4 , Cs 2 MnCl 4 , both being spin | 
compound, turn out to be a layered Heisenberg AFM 
with Ising like anisotropy [SS] . In this case the BKT sce¬ 
nario does not hold. On the other hand applicability of 
the Berezinskii-Kosterlitz-Thouless scenario on triangu¬ 
lar chromium-lattice AFM with S' = 3/2 has been in¬ 
vestigated only via electron spin resonance (ESR) tech¬ 
nique |69| . In this regard it is worthwhile to mention 
that the spin dynamics in a based spin | Hon¬ 

eycomb lattice anti-ferromagnetic material MnPSs has 
been investigated via INS experiment m- The critical 
properties of this material have been reported to be well 
described by 2D XXZ Hamiltonian (anisotropy param¬ 
eter being equal to 0.998) only in the low q* regime. 
We have performed an initial study towards the calcula¬ 
tion of DSF for MnPSs using our model and this reveals 
that the semi-classical BKT phenomenology is produc¬ 
ing a much better agreement with the experimental data 
corresponding to this material at 85 K m- The agree¬ 
ment is only in terms of peak position and the range of 
energy transfer over which we get acceptable (positive) 
values of DSF. This range encompasses a considerable 
portion of the range of experimental interest. Regardless 
of this, the possibility of meron dynamics in this material 
can not be adequately described with the available INS 
data. To be very specific, the available INS data is at 
85 K whereas, the critical dynamics is completely con¬ 
fined within the plane (and therefore, corresponding to a 
perfect 2D system) only above 105 K (TOj. Furthermore, 
in this material Tpf > Tbkt and therefore, above Tbkt 
and below T^v well defined magnon modes persist which 
is in contrast to the case corresponding to La 2 Cu 04 . 

Presently it seems that there is a scarcity of INS 
data in search for the spin dynamics induced by topolog¬ 
ical excitations corresponding to XY anisotropic (easy 
plane anisotropic, i.e., XXZ type) layered Heisenberg 
anti-ferromagnetic materials having higher values of spin. 
More INS experiments on this type of materials (in a 
suitable temperature range) would be quite interesting 
in view of the fact that it would serve as a good testing 
ground for the applicability of such a conventional semi- 
classical theory of spin dynamics induced by topological 
excitations. 

It is also very important to emphasize the fact that 
in our calculational analysis the structure of the classical 
vortex/meron has been built in the background of the 
Nfel state (see (§ of section H). Since the Neel state is 
not an exact ground state for the two-dimensional quan¬ 
tum anti-ferromagnetic spin systems, such a choice fur¬ 
ther adds to the reasons for the above mentioned unphys¬ 
ical behaviour. 

Keeping aside the occurrences of unphysical nega¬ 
tive values of the integrated intensities obtained from the 
conventional semi-classical-like theory, the results agree 
with those from the experiment quite well in terms of 
the existence of the central peak. However, this peak is 


a genuine one only in the sense that it occurs well within 
the experimental resolution width. 

Moreover, quantitative disagreement, at the finite 
values of energy transfer, between the magnitudes of the 
theoretically calculated integrated intensities and those 
obtained experimentally strongly indicates that forma¬ 
tion of the topological excitations from fragile collective 
modes and interactions between them can play a very 
crucial role in the spin dynamics. The values of the mo¬ 
ment at different temperatures further strongly suggest 
the same. Therefore, the interplay between the frag¬ 
ile single magnons or multi-magnon composites and the 
topological excitation is very crucial for the proper un¬ 
derstanding of the dynamics induced by free movement 
of the latter in the quasi-two-dimensional low-spin anti¬ 
ferromagnetic systems. This further takes care of the 
fragile magnon contributions to the total integrated in¬ 
tensities as well which has not been filtered out from 
the experimental data (for the reasons stated earlier in 
the section HI). In addition, the severe inadequacy of 
the semi-classical treatment in the case of quantum anti- 
ferromagnet may also contribute to such quantitative dis¬ 
agreement for the integrated intensities. 

Our investigations presented in this communication 
establish the fact that a complete quantum treatment 
is essential to describe the detailed features of the dy¬ 
namics of mobile topological excitations corresponding 
to the quasi-two-dimensional low-spin anti-ferromagnetic 
systems, by taking into consideration the vortex/meron- 
fragile magnon interactions as well. This further brings 
out the possibility of a characterization of different lay¬ 
ered materials into two classes, viz., conventional BKT- 
like and quantum corrected BKT-like. Moreover, since 
the Neel state is not an exact ground state for the two- 
dimensional quantum spin systems, the construction of 
the quantum vortices/merons in the complete quantum 
treatment has to be supplemented with proper choice of 
the ground state too. However, calculation of the dynam¬ 
ical structure function in a completely quantum mechan¬ 
ical framework is highly non-trivial. Earlier attempts 
towards this goal could not explain the occurrence of the 
“central peak” in the DSE obtained in the INS exper¬ 
iment performed on several quasi-two-dimensional ma¬ 
terials [SI [72H71]. Quite recently a theoretical frame¬ 
work had been developed based on the spin coherent 
state path integral formalism to describe the topolog¬ 
ical properties of static vortices and anti-vortices |7S1- 
l80| . An extension of this framework to the case of mo¬ 
bile spin vortices/merons (and anti-vortices/anti-merons) 
is crucial for quantum mechanical calculations of the 
DSE and the integrated intensity as well. Insights from 
Quantum Monte Carlo calculations may be quite use¬ 
ful in this endeavour m- This would be taken up 
in future. More importantly the short range 2D anti¬ 
ferromagnetic spin-spin correlation persists even in the 
superconducting phase of the underdoped cuprates. This 
has been investigated and verified via INS experiments 
[82] . Moreover it has been shown very recently that the 
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spin fluctuations in the AFM quantum critical region of 
the Fe based superconductors and some heavy fermion 
compounds can be modelled by dissipative quantum XY 
model and hence, the static and dynamics of topological 
excitations are the key factors for 2D spin-spin correla¬ 
tions [SH [85]. Therefore, this investigation of the dy¬ 
namics of the BKT vortices/merons may also contribute 
substantially to the microscopic understanding of lightly 
doped anti-ferromagnetic cuprates and the above men¬ 
tioned other systems as well [471 IS3H85] . 
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Appendix A: The Tukey and the modified Tukey 
function 

The most general form for the Tukey function is given 

by, 


R{t) 


1 TT 2t TT 

-[l-kcos(--- -k --- -tt) 

^ i Oi tvyi -L CX- 


1 . , TT 2t TT , , 

x[l + cos(—— -7— 

^ i Qi Cttt, -L Ck 

vO , 


for - < t < - a 

2 - 2 

r tm ^ ^ 

2 - - 2 

r ^ 

for a < t < — 

2 “2 

otherwise. 


(Al) 


where a is called the tapering parameter [sniini]- 

The Tukey function we have used in this communi¬ 
cation is corresponding to a = 0. In this case the above 
general expression takes the form. The Tukey function is 
given by. 


R{t) 


1. /SttF, 

-[1-f cos( —)], 

. 0 , 


for |f| < 
otherwise. 


(A2) 


which is corresponding to zero tapering [601 161] . The 
Fourier transform of the above Tukey function (corre¬ 
sponding to zero tapering) is given by. 


, 1 . / ^t-m , / 2 1 

R{u}) = -^ 




). (A3) 


The full width at half maximum (FWHM) for the above 

/T-'\ . 

function is given by expressed in the 

units of energy, where the superscript T signifies the 


Tukey function [53]. To find the value of tm the above 
expression for is equated to the value of the 

resolution width, 1.4 meV at FWHM, specified in the 
experimental observations corresponding to La^CuO^. 

The modified Tukey (MT) function we have used in 
our analysis is corresponding to 50% tapering and it is 
given by [BO] El] , 


nit) = < 


r , 47rt , , tvji 

[1 — cos(-)], for —— < t < 


t<n 


r tm , , , 

for - < t < - , . , 

4 4 (A4) 

tm. . . . 


i[l-cos(^)], forff<t<^ 


2 

10 , 


otherwise. 


Fourier transform of the above modified Tukey function 
is given by. 


)K- - 




(A5) 
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The full width at half maximum corresponding to the 
above function 'R.{uj) can be found out to be « 


, expressed in the units of energy. In this case, to find 

the value of tm we have equated to the value of 

the resolution width specified in the experiment. 
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